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EFFECTS  OF  MIRROR  COATING-THICKNESS  NONUNIFORMITIES 
ON  POLARIZATION  AND  OPTICAL-FIGURE  PROPERTIES  OF  LASER 
AND  BEACON  LIGHT 

Xiong  Shengraing  and  Zhang  Yundong 

Institute  of  Optics  and  Electronics 
Chinese  Academy  of  Sciences 
P.O.  Box  350 

Shuangliu,  Chengdu  610209 


ABSTRACT  This  paper  discusses  mirror  coating  nonuniformitics  caused  by  geometrical  con 
straints  in  the  chamber  used  and  their  effects  on  polarization  of  light  beams  for  enhanced  metal 
coa ti ngs.  Phase  distortions  of  1.06  fim  laser  beam  and  0.5893  —  /im  beacon  beam  have  been  calcu 
lated.  Phase  distortions  with  same  signs  have  beer  abtaincd  using  a  modined  design. 

KRY  WORDS  coaling  thickness  nonuniformitics  ,  polarization,  phase  distortion. 


I .  Introduction 

The  mirror  diameter  is  D  at  the  beginning  teimiinal  of  a  high- 
power  laser.  Based  on  considerations  of  Gaussian  light  beams,  at 
the  beam  waist  the  divergence  angle  6  is  approximately  equal  to  A/D 
(A  is  the  laser  wavelength).  Also,  consideration  is  given  to 
thermal  blooming  generated  by  a  high-power  laser  in  atmospheric 
propagation,  and  to  other  nonlinear  optical  effects,  such  as 
stimulated  Raman  scattering;  it  is  required  that  the  mirror 
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diameter  D  should  be  large  at  the  beginning  terminal .  In  a 
practical  high-power  laser  emission  system,  the  mirror  diameter  is 
generally  between  5  and  20m.  A  large-diameter  reflective  mirror  is 
not  only  difficult  to  machine,  but  is  also  high  in  cost. 
Additionally,  depositing  a  thin  coating  on  this  large  mirror 
surface,  especially  the  uniformity  of  coating  layer  thickness 
distribution,  is  very  important  and  difficult.  Nonuniform 
deposition  leads  to  polarization  aberrations  and  wavefront 
anomalies,  thus  affecting  far  field  diffraction  properties  [1]. 

Another  important  problem  in  a  high-power  laser  system  is  the 
presence  of  many  factors  leading  to  optical  wavefront  errors  in  the 
optical  passage  so  that  the  light  beam  energy  is  scattered.  Thus, 
compensations  should  be  made  with  adaptive-optics  technology. 
However,  an  appropriate  light  source  should  be  the  information 
carrier  of  the  wavefront  errors  for  light  from  the  source  in  order 
to  compensate  for  these  errors .  Such  light  is  usually  called  a 
beacon  light. 

After  the  light  beam  bounces  back  from  the  reflective  mirrors 
of  different  coating  systems,  phase  shifts  will  be  generated. 
Light  at  different  wavelengths  undergo  different  phase  shifts. 
However,  due  to  the  nonuniformity  of  coating  thickness,  phase 
anomalies  are  generated.  With  respect  to  light  at  different 
wavelengths,  phase  anomalies  are  different,  even  generating  phase 
anomalies  with  sign  reversal.  Then,  compensating  with  a  beacon 
light  may  yield  reverse  results,  failing  to  attain  the  expected 
purposes . 
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From  the  available  conditions,  this  article  takes  into  account 
the  geometric  constraints  in  the  film-coating  chainber  in 
calculating  the  coating-thickness  distribution  of  a  nearly- 
spherical  surface  reflective  mirror  with  0.8m  diameter  in  order  to 
attain  optimization.  In  the  case  of  the  metal-enhanced  type  high- 
reflectivity  mirror  system  of  this  distribution,  with  respect  to 
the  phase  anomalies  generated  by  a  1.06|jim  primary  laser  and  a 
0.5893|im  beacon  light,  and  with  the  modified  coating  system,  the 
design  is  so  executed  that  the  phase  anomaly  signs  at  both 
wavelengths  are  consistent. 

II.  Effect  of  Physical  Constraints  in  Film-coating  Device  on 
Coating-Thickness  Uniformity 

Two  kinds  of  uniformity  should  be  of  concern  to  us  for 
coatings  on  large  reflective  mirrors:  one  is  orientational 
uniformity  and  the  other  is  diametral  uniformity.  Orientational 
uniformity  can  be  very  easily  attained  by  high-speed  rotation  of 
the  reflective  mirror,  or  compensated  with  a  ring  vapor  source. 
Fig.  1  shows  the  ring  vapor  source.  With  regard  to  large 
reflective  mirrors,  for  safety  reasons  generally  low-speed  rotation 
was  employed.  Therefore,  it  was  necessary  to  use  a  ring  source 
apparatus  for  good  orientational  uniformity.  However,  diametral 
coating-thickness  distribution  was  difficult  to  control.  Diametral 
coating-thickness  distribution  relied  on  three  parameters [ 2 ] : 

(1)  the  f-value  of  an  optical  reflective  mirror  and  surface 
aspherical  shape;  (2)  vapor  source-to-substrate  distance  H; 
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Fic  .1  Riny  source  cuahnp  geometr)  using  a  sen  e‘> 
of  vapor  source  arranged  lo  simulate  a  ring 
vapor  source  of  diameter  2D  placed  endern- 
caib  a  rotating  primary  mirror. 

1.  rotating  primary  mirror 

2.  mirror  diameter  D, 

3.  vapor  stream 

4. sourcc  -  to  -  subtraic  distance  H 

5.  source  ofTsei  D  6.  vapor  source 


Radial  distance  (m) 

Fig -2  ( -> 01 'm coaling  ihicKncv^  unirormit\  fo' 

O.S  meter  diameter  F' 2.0  concave  mtrror. 
i  .mirror  ;;  epomum  source’  “  lo-  suntraic  hcigr.: 

//  =  0.9m  3. optimum  ofTscl  £>=  0.63m 

4  \apo:  source 


3.nci  (3)  distanc©  D  (offs©t)  of  th©  vapor  soiiirc©  ir©lativ©  to  th© 
rotational  center  of  the  reflective  mirror.  With  respect  to  the 
spherical-surface  rotation  fixture,  the  coating  thickness 
relationship  at  any  point  P  on  the  substrate  is: 


_  mh 

[  (/i*  -i-  Z>^  -r  )  (/icos/i  T  psinA  )  —  2 p  lysinA] 

( 1 ) 

*  i  p 

In  the  equation 

/)—//  —  (  1  *“  COSyd  ) 

(2) 

p 

sm^=  — 

(3) 

In  the  above  equations,  m  is  the  total  quantity  of  evaporated 
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material;  p.  is  density;  R  is  the  spherical-surface  radius;  p  is  the 
ratio  between  the  t^p  (relative  thickness)  and  thickness  in 
Eq.  (1)  when  p=0,  from  point  P  to  the  diametral  radius  of  the  axis 
of  rotation. 

Fig.  2  indicates  the  optimal  coating-thickness  condition  for 
a  reflective  mirror  with  F/2.0,  and  an  0 . 8m-diameter .  From  this 
figure  we  can  see  that  the  height  of  the  film-coating  apparatus 
housing  had  to  be  1.3m  and  the  diameter  had  to  be  1.5m  for 
uniformity  below  0.5%  with  the  vapor  source-to-substrate  distance 
H  0 . 9m  and  the  offset  0.3m.  However,  the  equipment  available  to  us 
could  not  satisfy  these  requirements.  Due  to  reasons  of  cost, 
remodeling  of  the  available  equipment  was  necessary;  thus,  ideal 
coating-thickness  uniformity  could  not  be  realized.  Fig.  3 
indicates  the  effect  of  offset  D  on  coating-thickness  uniformity 
with  known  height  H.  From  the  figure,  for  H  at  0.65m,  for 
uniformity  of  ±1.5%,  the  offset  was  required  to  be  between  0.48m 
and  0.5m.  However,  our  available  coating  machine  was  1.5m  in 
diameter;  a  maximum  offset  of  0.46m  could  be  attained.  Fig.  4 
indicates  the  0.45m  offset.  With  respect  to  coating-thickness 
variation  of  high-order  functions,  the  figure  indicates  that 
uniformity  in  height  ±2%  can  be  obtained  between  0.56  and  0.64m. 

In  the  foregoing  discussion,  the  cover  technique  was  not 
applied.  For  unifoirmity  calculations,  it  was  assumed  that  all 
vapor  sources  were  small-area  sources ,  obeying  the  cosine 
distribution  law.  Spherical-surface  treatment  was  carried  out  on 
the  reflective  mirror. 
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Fig.  3  Coaiing  thickness  non  -  uniformity  for  O.Smet  cr 
diameter  F/2.0  concave  mirror  for  various  offset 

distance  and  for  a  fixed  source— to  “  substrate  di¬ 
stance  of  0.65  meters. 


Fig  .4  Coating  thickness  non- uniformity  for  0.i< 
meter  diameter  F/ 2.0  concave  mirror  disf 
ance  for  various  height  and  a  fixed  offset 
distance  of  0.45  meters. 


III.  Initial  Design  of  Coating  Thickness 

The  authors  selected  the  wavelength  of  the  primary  laser 
wavelength  in  the  neighborhood  of  l|im;  l.OSjim  was  decided  on.  The 
beacon  light  could  be  the  light  diffusely  reflected  from  the  target 
by  the  primary  laser,  or  a  laser  at  other  wavelengths  in  the  relay 
satellite.  The  authors  employed  an  artificial  beacon  light 
generated  by  resonance  scattering  of  an  0.5893M.m  laser  by  sodium 
atoms  in  the  alkali-metal  layer  in  meteors  in  the  80-100km  upper 
atmosphere,  as  proposed  by  a  French  astronomer,  Labyrie.  It  was 
assumed  that  reflectivity  greater  than  99%  was  required  for  1.06|im 
wavelength;  however,  reflectivity  greater  than  95%  was  required  for 
0.5893|im  wavelength. 
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The  initial  coating  system  design  was 

G/Cr/Ag/ AJ,03/(LH)VAir 

L  and  H  represent,  respectively,  the  Si02  and  the  TazOs  coating  at 
one-quarter  wavelength  optical  thickness  in  the  primary  laser 
wavelength.  Cr  and  AI2O3  were  the  adherent  layers;  G  was  the 
substrate;  and  the  geometric  thickness  of  Ag  was  70nm. 

A  design  was  selected  that  deposited  the  full  one-quarter  or 
ha  If -wave  length  medivim  coating  onto  the  metal;  this  was  chosen 
because  of  easy  monitoring.  The  (LH)^  structure  was  applied  in  the 
initial  design;  this  was  so  because  of  wider  reflectivity  at  both 
wavelengths  of  interest  to  the  researchers .  Ag  was  selected 
because  it  has  the  high-reflectivity  L-spectral  band,  from  visible 
light  to  infrared.  From  this  design,  we  can  see  that  the  maximum 
incident  angle  of  the  primary  reflective  mirror  was  15°.  In  the 
following,  a  discussion  is  presented  on  the  polarization  properties 
(ellipticity  and  phase  anomalies)  due  to  nonuniform  coating 
thicknesses  with  respect  to  the  two  wavelengths  of  interest. 

IV.  Effect  of  Thickness  Nonuniformity  on  Ellipticity  in  the  Initial 
Coating  System  Design 

When  linearly  polarized  light  is  incident  on  the  primary 
reflective  mirror,  other  than  in  the  normal  direction,  the  combined 
reflected  light  is  elliptically  polarized  light.  In  other  words, 
there  should  be  a  phase  difference  in  the  reflected  light  between 
the  light  vector  P  parallel  to  the  incident  surface,  and  the 
optical  vector  S  noimial  to  the  incident  surface.  Moreover,  there 
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were  different  reflective  oscillation  amplitudes  between  the  two 
mutually  orthogonal  components  of  the  two  phases.  Fig.  5  indicates 
a  functional  diagram  of  coating  thickness 
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Fig.  5  Maximum  di/Terence  in  p  and  j 
components  of  reflectivity  and 
phase  for  the  initial  selected 
design  as  a  function  of  thickness 
non -uniformity  for  15'  incidence 
angle . 


Rtiativc  non  —  uniformity  {%) 


nonuniformity  due  to  variations  of  refractivity  and  phase  of  the 
primary  laser  wave  and  the  beacon  light  wave.  At  1.06)im,  the  phase 
^p~As<l°.  Even  when  the  thickness  nonuniformity  was  4% 
and  the  reflectivity  difference  in  intensity  was  smaller  than 
0.02%,  the  difference  in  reflectivity  was  less  than  0.02%  at 
0.5893jim,  and  the  phase  difference  less  than  3.5°.  Generally 
speaking,  there  was  little  effect  on  ellipticity  from  the  4% 
thickness  nonuniformity  for  very  oblately  polarized  light  as  the 
reflected  light  with  a  15°  angle  of  incidence. 

V.  Effect  of  Coating  Thickness  Nonuniformity  on  Phase  Anomalies  at 
Reflected  Waves  at  Different  Wavelengths 
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Nonuniformity  in  thin-coating  properties  was  exhibited  as 
different  optical  thicknesses  and  phase  shift  differences  (phase 
distortions).  The  phase  shift  differences  will  lead  to  optical- 
figure  distortions,  especially  greatly  affecting  the  far-field 
diffraction  properties  of  the  system.  Fig.  6  indicates  the  phase 
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Fig  .  b  Phase  dislorlion  in  degrees  for  the  inilial 
scJccicc  de.sicn  ai  the  iwo  waveiengihs  a.s 
0  function  of  thickness  non  *“  uniformiiy 
for  O.S  meter  diameter  /*/2.0. 


Fig  .  7  Phase  distortion  for  the  initial  selected 
design  along  ine  radial  distance  of  0.i^- 
meter  diameter  F!  2.0  mirror  for  source 
“  io“substraic  height  of  0.55  meters 
and  offset  distance  of 0.45  meters. 


distortions  in  the  two  reflected  waves  at  the  wavelength  due  to 
coating-thickness  nonuniformity  in  the  initial  coating  system 
design.  From  the  figure  we  can  see  a  -11°  or  -A/32.7  wavefront 
distortion  can  be  generated  by  +4%  nonuniformity  at  1.06M.m;  a  phase 
shift  of  -2.4°  or  a  wavefront  distortion  of  A/150  are  generated  by 
1%  nonuniformity.  The  phase  distortions  were  greater  at  0.5893)im. 
The  phase  shift  of  140.28°  or  a  wavefront  distortion  of  A/257  was 
generated  by  4%  nonuniformity.  A  phase  shift  of  -8.5°  or  wavefront 
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distortion  of  A./42.3  was  genoxated  with.  1%  nonunif onnity . 

With  respect  to  a  reflective  mirror  of  F/2.0  and  0.8m 
diameter.  Fig.  7  shows  the  phase  distortions  of  diametral 
distribution  as  indicated  by  wavelength.  In  the  computation,  it 
was  assumed  that  the  large  mirror  surface  is  approximately 
spherical;  the  vapor  source-to-substrate  distance  was  0.55m;  and 
the  offset  vapor  source  was  0.45m. 

Sign  reversal  of  phase  distortions  will  lead  to  other 
problems.  Optical-figure  errors  at  a  given  wavelength  can  be 
partially  revised  with  adaptive-optics  compensation.  If  there  is 
sign  reversal  of  wavefront  distortions  of  a  light  beam  at  a 
wavelength,  then  the  wavelength  distortions  at  the  wavelength  may 
further  degrade  the  wavefronts  at  another  wavelength  during 
modification.  Wavefronts  at  the  two  wavelengths  should  be  modified 
at  the  same  time.  The  light  beams  at  these  two  wavelengths  applied 
the  same  phase  distortion  signs .  Phase  distortions  of  the  light 
beams  at  both  wavelengths  should  be  consistent,  if  possible. 
Therefore,  another  modification  design  was  considered. 

VI.  Effect  on  Polarization  Properties  and  Phase  distortions  in  a 
Modified  Design 

For  the  same  phase  distortion  signs,  numerous  simple  designs 
were  studied.  Based  on  theoretical  calculations,  the  modified 
coating  system  design  was 

G/AJ.Oj/Ag/Al.Oj/  (LH)~/ Air 

H  and  L  indicate,  respectively,  the  Ta203  coating  and  the  SiOz 
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coating  of  one-quarter  wavelength  thickness  at  the  modified 
wavelength  1.17}am;  AI2O3  is  the  adherent  layer.  Fig.  8  shows 


S  Maximum  difl'ercnce  in  p  and  s  components 
of  reflectivity  and  phase  for  modified  design 
as  a  function  ot  thickness  non  “  uniformity 
for  15'’  incidence  angle  . 


Fig  .  9  Phase  distortion  in  degrees  for  the  modified 
design  at  the  two  wavelengths  as  a  function 
of  thickness  non  -  uniformity  for  0  8  -  meter 
diameter  /'/ 2.0  mirror. 


the  effect  of  thin-coating  polarization  properties  due  to 
nonuniformity  in  design  thicknesses .  It  was  discovered  that  this 
thin  coating  system  design  had  greater  variation  with  thickness 
nonuniformity  in  reflectivity  at  0.5893|im  in  this  coating  system 
design.  Both  designs  had  corresponding  reflectivities  at  1.07|am. 
By  comparison  with  the  initial  design,  the  reflectivities  at  the 
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nonuni formity  value  of  +4%  were  slightly  higher  in  the  modified 
design;  the  reflectivity  in  the  modified  design  at  0.5893nm  was 
lower  than  in  the  initial  design.  In  the  modified  design.  Fig.  9 
shows  the  phase  distortions  generated  by  two  wave  beams  at 
different  wavelengths  due  to  different  nonuniformity.  From  the 
figure,  we  can  see  that  a  phase  shift  of  -23°  at  1.06)im,  or 
wavefront  distortions  of  A/15.6  were  generated  by  4%  nonuniformity 
at  1.06)im;  a  phase  shift  of  -4.5°  or  wavefront  distortions  of  A/80 
were  generated  by  1%  nonuniformity.  Distortions  were  greater  at 
0.5893)im.  A  phase  shift  of  -55°  or  wavefront  distortion  of  A/6.5 
were  generated  at  4%  nonuniformity.  A  phase  shift  of  -14.5°  or 
wavefront  distortions  of  A/24.8  were  generated  by  1%  nonuniformity. 
For  nonuniformities  at  different  thicknesses,  there  were  the  same 
signs  of  light  beam  wavefront  distortions  at  the  two  wavelengths. 


0.02 j  Radial  distaocc  (m) 


Fig.  10  Phase  distortion  for  the  modified  design 
along  the  radial  distance  of  0.8 -meters 
diameter  f/ 2.0  mirror  for  source- to- 
substrate  height  0.56  meter  s  and  ofTsc: 
distance  of  0.45  meters  . 


Fig.  10  presents  a  reflective  mirror  of  F/2.0  and  diameter 
0.8m;  in  the  modified  coating  system  design,  there  are  phase 
distortions  of  the  two  wave  beams  along  the  diametral  distribution. 
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VII.  Conclusions 


From  the  foregoing  analysis,  in  the  modified  coating  system 
design  (by  depositing  optical-thickness-medium  coatings  of  one- 
quarter  or  one-eighth  wavelengths  onto  metallic  silver),  there  are 
the  same  phase  distortion  signs  generated  by  light  beams  at  the  two 
wavelengths  due  to  coating  thickness  nonuniformities.  At  15° 
incidence,  there  is  a  smaller  effect  on  polarization  in  the 
modified  design.  It  is  of  interest  that  the  reflectivity  was  not 
reduced  in  the  wavelength  range.  This  design  was  suitable  for 
large  mirror  surfaces  0.8m  in  diameter.  For  easy  monitoring,  all 
medixim  coatings  were  controlled  at  the  same  wavelength.  Such  same- 
wavelength  monitoring  reduced  the  property  changes  due  to  chromatic 
dispersion  of  each  coating  layer  under  multiple-wavelength  control. 

The  first  draft  of  the  article  was  received  on  July  23,  1991; 
the  final,  revised  draft  was  received  for  publication  on  September 
20,  1991. 
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EFFECT  OF  PHASING  ERRORS  ON  BEAM  COMBINATION 


Hu  Zhiping,  Le  Shixiao,  Song  Ruhua,  and 
Zhou  Hongyang 

Hu,  Le  and  Song  of  University  of  Electronic  Science  and  Technology, 
Chengdu  610054;  Zhou  of  China  Academy  of  Engineering  Physics,  P.O. 
Box  532,  Chengdu  610003 

BSTRACT  Based  on  the  research  beam  control  for  synthetic  aperture  systems  in  free  space. 

■Jis  paper  analyzes  and  discusses  efTecis  of  phasing  errors  on  beams  combination.  The 
basing  errors  include  piston  errors,  till  errors  and  focus  errors.  The  corresponding  physical 
lodels  arc  developed  .Finally  results  of  numerical  computations  arc  given. 
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I .  Introduction 

For  a  synthetic-aperture  system  as  a  single  coherent  high- 
power  laser  emission  system,  in  recent  years  the  system  was  of 
growing  concern  to  researchers .  By  controlling  the  wavefront 
phases  of  beams  (emitted  at  various  subapertures)  of  a  synthetic- 
aperture  system,  coherent  far  field  combination  of  light  beams  can 
be  carried  out.  However,  when  controlling  the  wavefront  phases  of 
the  emitted  light  beams,  generally  there  will  be  errors  in  phase 
control  with  various  factors  affecting  phase  control,  leading  to 
degraded  properties  of  the  synthetic -aperture  system  and  reducing 


the  peak-value  light  intensity  coherently  combined  in  the  far 
field.  The  authors  developed  studies  on  the  concept  of  light  beam 
control  of  the  synthetic-aperture  system  in  the  ideal  situation 
without  considering  various  factors  affecting  precision  phase 
control  of  the  light  beams.  Thus,  a  physical  model  for  coherent 
and  noncoherent  combination  of  light  beams,  as  well  as  beam¬ 
pointing  and  beam-shifting  were  are  presented.  Theoretical 
analysis  and  niimerical  simulation  calculations  were  conducted 
[2,3].  On  the  basis  of  the  developed  work,  the  article  further 
analyzes  and  discusses  the  factors  affecting  precise  control  of 
light  beam  phases  (piston  and  tilt  errors  of  tuned  phase,  and 
tuned-focus  errors),  and  the  effects  on  coherent  combination  of 
light  beams.  In  addition,  numerical  simulation  calculations  and 
three-dimensional  diagrammatic  demonstrations  were  conducted. 

II.  Theoretical  Analysis  and  Numerical  Simulation  Calculations 

For  the  theoretical  analysis,  the  authors  assumed  that  the 
light  beams  emitted  from  the  emission  system  of  a  laser  device  are 
plane-wave  laser  beams ; 

V, )  (x,  ,  >», ) 

In  the  equation,  B  is  a  constant.  For  convenience  in  the  following 
analysis,  we  set  B=l;  pi(Xi,  y^)  is  the  pupil  function  of  the 
emission  system,  with  the  following  definition: 

0  (.tj+v?)  >(rf/2)= 
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In  the  equations,  d  is  the  pupillary  diameter. 


Assume  that  the  synthetic-aperture  system  is  composed  of  N=lxm 
subapertures.  For  each  subaperture,  its  center  coordinates  are 
Yt(i/j)/i=lr2,...,l;  j=l,2,.. ..,m.  The  combined  light  field 
at  the  far  field  target  is[2]: 
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In  the  equations,  A(x,y,  z  )=exp[  (  jkz )  •exp(  (  jk/22: )  (x^+y^)  ]  /  (  jAz ) , 

F{  }  represents  the  Fourier  transform.  k=2TT/A,  X  is  wavelength, 
and  (Xt,  yt)  are  the  coordinates  of  the  pointing  dot.  Let  Xt=yt=0 
(the  light  beams  emitted  from  various  subapertures  converge  onto 
the  axial  line),  then  we  obtain 
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Obviously,  from  Eq.  (3)  we  know  that  the  peak-value  light 
intensity  in  coherent  combination  is 

/*,=  IC’(0,0  l3=-{p,(r,  )}P-(/xm)  =  =A^^  •  /,m  (4) 

In  the  equation,  I^^  indicates  the  peak-value  light  intensity 


of  a  single  subaperture  system  at  (P,0,z). 


In  the  following  analysis  and  calculations,  we  let  N=16,  d=lin. 


X=l|im,  and  z=10^kin. 
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Fig.  1  shows  the  light-intensity  distribution  of  beams 
coherently  combined  in  the  synthetic-aperture  system;  the  data  in 
the  figure  indicate  the  relative  values  of  peak-light  intensities, 
which  are  the  same  as  later  on.  Fig.  2  shows  the  regional  energy 
distribution  curve  of  light  beams  coherently  combined  in  the  far 
field.  In  the  figure,  Rp  indicates  the  radius  of  the  Airy  disc. 
E(AIRY  DISC)/E(TOTAL)=0.868299E+00  indicates  that  the  percentage  of 
energy  in  the  first  Airy  disc  is  86.8299%  of  the  total  energy. 
Rd=0.61  X  X  X  z/b;  the  unit  is  0.61  x  X  z/b.  X  is  wavelength;  z  is 
the  mirror-to-target  distance;  and  b  is  the  equivalent  pupillary 
radius  of  the  phase  control  array. 

1.  Effect  of  piston  errors  [4]  of  the  tuned  phase 


n 


The  piston  errors  of  the  tuned  phase  stem  from  the  phase-error 
term  when  tuning  and  controlling  the  wavefront  phase  of  the  light 
beams  emitted  from  the  subapertures,  such  as  exp(jkcj)ij)  (4)^^^  is  not 
related  to  Therefore,  the  light  field  distribution 

coherently  combined  in  the  far  field  varies  as: 
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Thus ,  the  coherently  combined  peak-value  light  intensities  become 
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In  other  words ,  the  piston  errors  of  the  tuned— phase  light  beams 
will  lead  to  a  decrease  in  peak-value  light  intensity. 

Figs.  3  and  4  are,  respectively,  the  light— intensity 
distribution  diagram  and  the  regional  energy  distribution  curve  of 
the  far-field  combination  light  beams  under  the  situation  of  (l)ij=0 
with  respect  to  4)jl=4>i2=4)22=0 . 5A .  From  these  figures,  we  can  see 
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that  the  piston  errors  in  the  tuned  phase  will  lead  to  decreased 
peak-value  light  intensities  with  light  beam  combination;  however, 
the  effect  is  smaller  on  energy  percentage  in  the  first  Airy  disc. 


2.  Effect  of  tilt  errors  of  tuned  phase  [4] 

Tilt  errors  of  the  tuned  phase  signify  the  phase-error  term 
due  to  tilt  and  wobbling  of  the  subapertures  in  bringing  forward 
such  as  exp{[  jk[Cijx(x-X3(i,  j)  )+Cijy(y-y3(i,  j)  )  ]}  when  controlling 
the  wavefront  phase  of  the  light  beams  emitted  from  the 
subapertures.  Here,  the  combined  light-field  equation  of  the  light 
beam  should  be 
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Obviously,  here  the  peak-value  light  intensity  becomes 
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Figs.  5  and  6  are,  respectively,  the  light-intensity 
distribution  diagram  and  the  regional  energy  distribution  curve  of 
light  beams  combined  when  Cijx=Cijy=0 . 5xl0'®rad.  From  the  above- 
mentioned  diagrams  and  curves  we  know  that  the  tilt  errors  of  the 
tuned  phase  not  only  reduce  the  peak-value  light  intensity  of  the 
far  field  beam  combination,  but  also  lead  to  energy  reduction  in 
the  first  Airy  disc.  The  greater  the  tilt  errors  in  the  tuned 
phase,  the  more  pronounced  is  such  effect. 


3.  Effect  of  light  beam  tuning  focus  error 

When  focussing  the  light  beams,  if  the  light  beams  are  not 
accurately  converged,  the  following  phase  term  of  focus  errors  [2] 
will  be  generated: 
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Fig.  7  is  the  light  intensity  distribution  diagram  at  z=5xl0'^m 
of  a  light-beam  coherently  combined  when  there  are  no  focus  errors. 
Fig.  8  is  the  light  intensity  distribution  diagram  while  the  beam 
combination  at  z=5xl0''m  when  there  are  focus  errors  (that  is,  when 

^11=^12=^21  =  ^22=^31=^32=^41=^42=0  •  1  ^^d  Ai3=Ai4=A23=A24=A33=A34=A43=A44  =  -0 . 1  ; 
in  the  physical  significance,  in  some  of  the  subapertures,  the 
light  beams  are  focussed  ahead  of  the  common  pointing  dot,  and  the 
remaining  subapertures  have  the  beams  focussed  behind  the  common 
pointing  dot.  Obviously,  we  know  that  the  peak-value  light 
intensities  here  are  apparently  reduced;  the  greater  the  focus 
errors,  the  more  pronounced  is  the  effect  on  the  beam  combination. 


III.  Conclusions 

In  the  foregoing,  detailed  analysis  and  discussion  were  given 
concerning  the  factors  affecting  precision  control  of  light  beam 
phases  in  the  synthetic— aperture  system;  moreover,  numerical 
simulation  calculations  were  conducted.  As  shown  in  the 
calculation  results,  since  precision  control  of  the  light 
beam  phases  cannot  be  carried  out,  this  will  lead  to  a  situation  in 
which  the  peak-value  light  intensities  of  the  beam  combination  are 


Fig.  7  The  intensity  distribution  of  coherent 
beams  combination  without  focus  errors 
(A,  =0).z  =  5x  10^  m. 

reduced.  However,  we  have  to  point  out  that  the  factors  affecting 
the  light-beam  phase  control  are  assumed  to  be  static  variables, 
not  random  and  dynamic  variables  in  the  above-mentioned  analysis 
and  discussion.  However,  with  respect  to  laser  propagation  in  the 
atmosphere,  of  more  practical  significance  is  a  discussion  of  the 
effect  of  beam-phase  precision  control  by  such  random  and  dynamic 
variables  of  atmospheric  turbulence.  In  later  work,  we  will 
gradually  develop  research  on  this  aspect  in  order  to  obtain 
physical  models  affecting  light  beam  phase  control  that  more 
practically  match  the  situation  in  atmospheric  turbulence. 
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ABSTRACT  Experiments  on  the  establishing  and  vanishing  time  features  of  thermal  lens  ,n  the 
mov.ng  or  static  absorbing  cell,  the  beam  diverging  features  at  far  field  through  the  blooming  cell 

the  comparision  of  the  beam  distortion  features  of  noruniform  absorption  with  that  of  uniform 
absorption  are  reported. 
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I .  Introduction 

When  a  high  power  laser  beam  propagates  through  the 
atmosphere,  air  along  the  beam  path  is  heated  since  the  gases 
absorb  laser  energy,  thus  inducing  the  thermal  blooming  effect  of 
the  laser  beam.  Research  on  the  thermal  blooming  effect  during 
the  atmospheric  propagation  of  high  power  laser  beams  was 
conducted  for  many  years  [1,2].  In  the  authors'  research 
problems,  a  laser  beam  will  be  propagated  upward  to  a  height  of 
about  1000km.  The  atmospheric  layers  are  thinner  compared  to 


this  1000-km  of  propagation  through  the  absorption  medium.  In 
addition  to  the  bulk  absorption  coefficient  directly  inducing  the 
thermal  blooming  effect,  there  is  also  the  effect  of  the 
atmospheric  wind  field  (including  light  beam  rotation).  After 
the  laser  beam  is  subjected  to  the  thermal  blooming  effect  when 
passing  through  the  absorptive  atmospheric  layers,  the  beam 
continues  to  propagate  for  a  long  distance  in  a  vacuum.  The 
variable  features  of  the  beam  are  the  subject  of  interest  to  the 
authors.  In  addition,  the  nonuniformity  of  atmospheric 
absorption  was  theoretically  analyzed  in  [3]  as  regards  the 
upward  propagation  thermal  blooming  effect  for  high  power  lasers. 
As  indicated  in  the  analytical  results,  there  are  obvious 
differences  between  uniform  absorption  and  nonuniform  absorption 
^^th  respect  to  the  thermal  blooming  effect.  In  other  words,  the 
^^^*^-^tional  concept  of  considering  the  atmosphere  as  consisting 
of  homogeneous  absorptive  layers  is  not  at  all  appropriate. 

Based  on  these  considerations,  the  article  reports  on 
simulation  experiments  that  the  authors  conducted  in  their 
laboratories :  ( 1 )  under  the  simulated  conditions  of  windy  and 
static  absorbing  gases ,  studies  were  conducted  on  the  time 
features  of  the  formation  and  disappearance  of  monopulse  laser 
thermal  lenses.  (2)  Under  conditions  of  uniform  and  nonuniform 
absorption  (the  total  absorption  coefficients  of  the  two  kinds 
are  kept  equal),  studies  were  conducted  on  the  differences  in 
thermal  blooming  of  a  laser  beam.  (3)  After  generation  of  a 
thermal  lens,  a  study  was  made  on  the  distortion  divergence 


properties  of  a  laser  beam  over  a  longer  distance  other  than  in 
an  absorptive  medium.  The  theoretical  calculations  in  the 
article  are  based  on  the  geometric-optics  method  described  in 
[1]. 

II.  Experimental  Conditions  and  Parameters 

Fig.  1  shows  the  entire  setup  of  the  experimental  system. 

The  primary  laser  was  a  TEACO2  monopulse  laser.  By  using  the 
P-Ge  pulling  device,  variations  of  the  central  light  intensity 
were  detected.  A  model  7623A100M  storage  [illegible] 
oscilloscope  was  used  in  recording.  He-Ne  light  and  a  CO2  laser 
beam  were  passed  coaxially  through  a  simulation  absorption 
channel;  a  CCD  detection  system  was  used  to  monitor  the  light 
beam  distortion  process.  A  four-section  structure  simulation 
absorption  channel  was  employed,  capable  of  passing  through  two 
light  beam  channels.  One  pass  (OD=30mm)  served  to  propagate  the 
high  power  laser  beam.  Another  pass  (OD=25mm)  was  used  to 
monitor  the  absorption  coefficient  in  the  various  sections  of  the 
channel  with  a  bomb  calorimeter  for  light  splitting  (into  faint 
light)  in  the  light  pass.  Propane  was  used  as  the  absorption 
gas .  This  gas  was  fed  while  controlling  the  amount  absorbed 
during  the  various  passes;  then  air  was  admitted  until  l.OlxlO^Pa 
was  reached.  In  the  figure,  the  concave  reflective  mirrors  Rj 
and  R2  tuned  the  10. 64m  laser  beam  into  a  nearly  aligned  light 
beam  to  pass  through  the  simulated  absorption  path.  The  light- 
spot  radius  at  the  1/e  power  level  was  5mm.  The  parallel  light 
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1.  synchronous  laser  power  2.  synchronous  signal  3.  simulated  absorption  system 
4.  gas  out  5.  gas  fill  6.  monitor  ‘7.  computer  sample  system  8.  PC  interface 
9  .  synchronous  control  10.  video  signal  I  l.CCD  camera  12.  view  screen 


tube  T  tuned  the  He-Ne  light  to  a  slightly  divergent  light  beam 
which  was  passed  through  a  germaniiim  lens  G  to  be  tuned  to  the 
same  axis  as  the  CO2  laser.  The  reflective  mirror  reflected 

the  He-Ne  light  several  times  in  the  laboratory  in  order  to 
obtain  the  intensity  distribution  after  passing  through  the 
medium  for  the  light  beam  in  far  field  propagation  (the  longest 
distance  was  z=30m) .  With  the  measurements,  the  authors  obtained 
the  following  experimental  parameters: 

(1)  TEACO2  energy  E=0.55J,  and  pulse  duration  t^=20^s 

(2)  overall  absorption  coefficient  a^^O. 193/m 

(3)  section-by-section  nonuniform  absorption  coefficient 

a  ,  =  0.4825 '  m,  a,  =  0.241 3/ m,  2.  =  0.0483 /rn,  2^  =  0. 

satisfying  {L/4)X2,  =  2,L 

(4)  length  of  absorption  pass,  L=4.6m 


(5)  thermal  distortion  parameter  N=  -In.  E^,  L-,  npc  a' =  \  es 
in  the  parameters  p=l .  2xl0'^g/cm^  (gas  density) 

Cp=lJ/g'K  (constant-pressure  specific  heat  of  the  gas 

•  5x10"^  K'^  (dn/dT,  derivative  of  refractivity  with 
respect  to  temperature) 

aco2=5mm  (light  spot  at  1/e  power  level) 

(6)  for  the  He-Ne  light  emitted  beam,  the  light-spot  radius 
ao=8.84mm  (at  channel  exit) 

a ( z=30 )=10 . 35mm  (30  mm  from  channel  exit) 

(7)  for  He-Ne  light,  the  Fresnel  nvimber  F( 30m)  =1.64 
(basically  satisfying  the  far  field. conditions ) . 

The  real-time  synchronous  quantized  CCD  detection  system 
consists  of  a  CCD-camera,  a  video  frequency  screen,  and  a 
computer.  In  the  central  portion,  the  main  indicators  of  the 
CAS 3  video  frequency  monitor  were  the  range  of  the  8 -bit  dynamic 
A/D,  and  D/A;  the  real-time  four-frame  512x512  memory  and 
25-f rames-per-second  acquisition  speed.  In  addition,  through  a 
74LS123,  two-channel  triggering  was  conducted,  thus  obtaining  the 
synchronous  triggering  of  the  video  frequency  field,  and  the 
time-division  delay  triggering  pulses  with  40ns  as  the  full  cycle 
in  order  to  carry  out  synchronization  between  the  CCD-camera  and 
the  monopulse  laser  device,  as  well  as  the  acquisition  of  light- 
spot  information  at  various  time  instants . 

III.  Data  Acquisition  and  Treatment  Method 

Due  to  the  restriction  on  the  dynamic  range  ( 8  bits )  of  the 


acquisition  board,  consideration  was  given  to  repeatability  of 
the  thermal  blooming  effect.  In  precision  synchronous  control, 
the  authors  applied  instant-by-instant  light  exposure  method  to 
obtain  the  nonsaturated  A/D  conversion  signal.  Fig.  2  shows  the 
one -dimensional  intensity  curve  of  the  original  He-Ne  light 


Measured  intensity  (9  bit) 


intensity.  By  using  two  light  exposure  times  of  l/250s  (pi  light 
intensity)  and  l/500s  ( p2  light  intensity),  the  authors  obtained; 

,Mp,-255  p,^255jA 

Po=  <  ‘ 

I  0  p.<255IA 


In  the  equation,  A=500/250=2,  pg  is  the  portion  of  saturated 
intensity  to  p^.  However,  the  actually  collected  light  intensity 
P“Pi+Po*  Thus,  the  measured  A/D  conversions  reached  the  maxijim 
value  of  512  (9  bits).  At  the  various  time  instants,  the  errors 
of  the  light-spot  intensity  e={ pds  does  not  exceed  5%. 
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Calculation  of  light-spot  radius: 


<r> 


ir  —  F  p  ds  ' 


r 

pds 

j 


(2) 


In  the  equation,  is  the  position  of  the  light-spot  center  of 
gravity;  ds=dxdy,  however 


ds=  dxdy ,  +  '  )  "'  *  M  x^  = 


xpds!  pds\y\=  ypdsj 


pds , 


With  respect  to  the  results  of  the  secondary  light  exposure, 
Eq.  (2)  can  be  processed  into: 


<r  >  = 


\r  -  r  p .  d  s 


\r-r^fp^d5 


p,as^ 


( 3  ) 


\p,ds 


Let 


p.ds-,E^=  p^ds:  <r:>  =^\r-r^fp,dsi  Ey,  <rl>  p^ds;  L^, 


We  obtain: 


<r2>  =  <r;>  +  <r\>EJ  E, 
\-^E,:E, 


(4  ) 


With  the  energy  statistics  of  the  light-spot  radii,  we  can 
calculate  that  the  light  beam  divergence  angle  is 

<fl>=A<r>;Ar 

In  the  equation,  A<r>  is  the  difference  between  two  points  of  the 
light-spot  radius  in  the  propagation  distance;  Az  is  the 
difference  of  the  distance  between  two  points. 
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IV.  Experimental  Results  and  Discussion 


Fig.  3  shows  the  variation  of  laser  beam  intensity  with  time 
at  the  center  and  after  distortion  of  the  COj  laser  beam.  From 
•  3  (b),  we  can  cloarly  see  the  obvious  nonlinear  attenuation 


(a)5/Js/div  (b)  5/is/div 

Original  intensity  Distort  intensity 

Fig  .  3  Curves  of  CO.  laser  intQnsity  at  centre  point 


effect  due  to  thermal  blooming.  In  other  words,  at  the  most 
crucial  instant,  t=7.5M.s  of  thermal  distortion,  the  light-spot 
center  intensity  is  nearly  zero.  This  results  from  expanding  the 
light  beam  with  the  thermal  lens.  Besides,  based  on  the  linear 
attenuation  difference  for  the  10.6)im  laser  due  to  the  absorption 
gas,  from  Fig.  3  we  can  see  that  the  time  t=5}i.s  when  the  thermal 
lens  begins  to  foirm. 

Fig.  4  shows  the  variation  of  the  thermal  blooming  light 
spots  at  various  instants  when  the  absorption  channel  is  quiet, 
and  after  the  action  of  the  CO2  laser  by  the  coaxial  He-Ne  light, 
due  to  the  thermal  lens  effect.  At  t=0,  this  indicates  the 
instant  of  CO2  laser  action.  We  can  see  that  in  the  early  period 
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After  pulse  f  =  40ms  f=80ms  Original 

(a)  Blooming  patterns  of  He-Ne  laser 
(a)  He-Ne 


Fig  .  4  Blooming  intensity  { 1 D  )  of  He  -  Ne  laser  at  diflerent  time 

of  laser  action  (t<20ms),  the  He-Ne  light  undergoes  a 
considerable  expansion.  The  light-spot  energy  radius  is  three 
times  the  radius  of  the  original  light  spot  at  the  1/e  power 
level.  This  result  is  quite  different  from  the  authors' 
theoretical  calculations  [5].  However,  in  the  later  period 
(t>40ms)  of  the  thermal  lens,  the  experimental  results  and  the 
theoretical  calculated  results  are  relatively  close  to  each  other 
(refer  to  Fig.  5).  This  indicates  that  the  later  period  of  the 
thermal  lens  has  the  thermal  propagation  playing  a  major  role. 
Under  such  conditions,  complete  restoration  of  the  light  spot 
requires  more  than  20ms. 

Fig.  6  shows  the  measured  He-Ne  light  when  the  absorption 


J3 


Normalized  intcnsily  2:  r- After  pulse  Nnrniaiizcti  »n‘ens>ly 


Fig  .5  Curves  calculated  on  blooming  intensii)  Fig  .  6  Time ~  variation  curves  measured  at  3-po'nis 

(ID)  of  He-  Ne  laser  condition  of  simulated  wind 

channel  moves  with  respect  to  the  time-based  variation  (near¬ 
field)  of  light  intensity  at  three  points  of  the  windward 
direction  radius.  Curve  A  is  the  light-spot  center;  B  is  the 
windward  point  r=a  beginning  from  the  center;  c  shows  the  reverse 
direction  r=-1.076a.  From  the  figure,  we  can  see  that  the 
restoration  time  of  point  A  is  90ms.  In  the  distortion  process 
of  point  B,  due  to  the  effect  of  the  wind,  the  light  intensity 
exceeds  the  original  value  of  light  intensity  at  the  center  at 
approximately  40ms.  The  restoration  time  is  150ms.  The 
foregoing  results  explain  that  the  atmospheric  wind  field  has  a 
certain  function  in  smoothing  out  the  thermal  lens  effect. 

However,  it  simultaneously  brings  with  it  the  problem  of 
deviation  of  the  light-beam  center  intensity. 

In  the  theoretical-calculation  model,  both  the  CO2  laser  and 
the  He-Ne  laser  applied  the  ideal  limited  gaussian  light  beams. 
However,  in  the  experiments  the  light-spot  intensity 
distributions  of  the  CO2  and  the  He-Ne  lasers  showed  certain 
fluctuations,  see  Fig.  2  and  [4].  With  respect  to  the  results  of 


the  experiments  and  the  theoretical  calculations,  there  are 
obvious  differences  when  t  is  less  than  tens  of  milliseconds. 

The  authors  analyzed  that  this  is  due  to  the  following  reasons 
[5]:  (1)  nonuniformity  of  the  CO2  intensity  distribution  leads  to 
fluctuations  in  refractivity  variation  in  the  propagation  path, 
thus  not  only  leading  to  enhancement  of  thermal  scattering  of 
light  intensity  in  the  small  scale  nonuniformity  of  the  CO2 
laser,  but  also  enhancing  the  nonuniform  light  intensity  in  the 
small  dimensions  of  the  He-Ne  light.  Thus,  the  scattered  energy 
to  a  greater  angle  is  greater  than  the  calculation  results  of 
using  the  homogeneous  gaussian  model.  In  other  words,  the 
observed  light  beam  expansion  is  greater  than  the  thermal 
blooming  expansion  of  the  light  beam  without  perturbation.  (2) 
Since  there  are  fluctuations  in  intensity  and  phase  for  the  He-Ne 
laser  beam,  under  the  action  of  the  coaxial  high  power  laser, 
light  beam  distortion  is  more  serious  than  distortion  of  the 
He-Ne  laser  beam  in  the  absence  of  perturbation  of  the  entire 
beam;  the  expansion  is  also  greater.  In  short,  the  experimental 
results  of  the  light  beam  expansion  are  apparently  greater  than 
the  theoretical  calculation  results,  for  the  following  reasons; 
small-dimension  light  intensity  and  phase  perturbation  bring 
about  instability  of  thermal  blooming.  However,  over  a  longer 
time  period  after  the  action  of  the  pulsed  laser,  the  thermal 
lens  tends  to  smooth  out,  and  the  effect  of  small-dimension 
nonuniformity  is  reduced.  The  results  of  the  experimental  and 
theoretical  calculations  of  light  beam  thermal  blooming  in  the 


absence  of  perturbation  tend  to  be  consistent. 


As  shown  in  Fig.  7,  after  the  action  of  the  monopulse  CO2 
laser,  the  light-spot  radius  of  the  near-field  distortion  of  the 
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Fig  .  7  Comparison  of  Blooming  radius  due  to 
nonuniform  absorption  10  that  due  to 
uniform  absorption 

□  nonuniform  absorption  uniform  absorption 


/  (ms) 

Fig  .  8  Time  -  variation  of  relative  sprea  ding—  angle 
of  He— Ne  beam  in  the  far  —  field 


coaxial  He-Ne  light  beam  varies  with  time.  The  upper  and  lower 
curves  indicate,  respectively,  the  results  of  nonuniform  and 
uniform  absorption  conditions.  In  the  experiments,  it  is 
maintained  that  the  overall  absorption  rate  ^a(z)d''=a  with  L 

remaining  constant  in  both  situations .  We  can  see  that  time 
periods  for  the  formation  and  disappearance  of  the  thermal  lens 
are  consistent  in  time.  However,  the  thermal  distortion  light- 
spot  radius  under  nonuniform  absorption  conditions  is  apparently 
larger  than  the  thermal— blooming  light-spot  radius  in  the  case  of 
uniform  absorption.  At  the  beginning,  the  difference  between  the 
two  cases  was  0.4mm.  In  other  words,  in  the  case  of  nonuniform 


absorption  conditions,  the  divergence  angle  of  the  thermal- 
blooming  light  beam  is  greater  than  the  divergence  angle  of 
thermal -blooming  light  beam  under  uniform  conditions,  by  about 
3%.  This  is  due  to  the  integration  relationship  of  thermal 
distortion  parameter  N  and  ^  'e.xp  [  -  a  ( .  Since 

the  propagation  direction  of  the  light  beam  is  from  high  to  weak 
absorption,  in  the  front  end  of  the  simulation  channel,  parameter 
N  is  greater  in  the  nonuniform  absorption  situation,  and  its  lens 
effect  is  also  greater.  However,  propagation  to  the  rear  end, 
since  the  reduced  absorption  of  laser  energy  and  the  reduction  in 
energy  distribution  density  due  to  distortion,  the  effect  of 
parameter  N  is  weakened.  Thus,  we  can  say,  upon  considering  the 
problem  of  upward  propagation  of  a  high  power  laser  through  the 
atmosphere,  it  is  inappropriate  to  consider  the  effect  of  thermal 
blooming  with  the  equivalent  uniform  absorption  model .  In 
addition,  due  to  the  reverse  directions  of  upward  and  downward 
propagation,  the  thermal  blooming  effect  is  irreversible  with 
respect  to  the  nonuniform  absorbing  atmosphere. 

Another  problem  we  examined  was  that,  relatively  speaking, 
the  atmospheric  layers  are  very  weak  for  upward  laser  propagation 
to  approximately  1000km.  With  respect  to  the  problem  of  very 
weak  absorption  layers,  how  the  generated  thermal  blooming 
affects  the  far  field  propagation  of  a  laser  beam  in  a  vacuum? 

Fig.  8  indicates  the  measurement  by  the  authors  under  the 
approximate  far  field  conditions  outside  of  30m  from  the  light 
channel  exit  in  the  experiments,  with  respect  to  the  results  of 
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the  time-dependent  ratio  between  the  thermal -blooming  light-beam 
divergence  angle  and  the  laboratory  light-beam  diffraction  angle. 
We  can  see  that  in  the  beginning  instant  of  thermal  blooming,  the 
light-beam  divergence  angle  due  to  the  thermal  blooming  effect  is 
nearly  eight  times  greater  than  the  diffraction  angle.  The 
result  is  basically  consistent  with  the  conclusion  when  the 
authors  calculated  under  the  conditions  of  the  corresponding 
parameter  N,  that  there  is  one  order  of  magnitude  in  the  ratio 
between  the  divergence  angle  and  the  diffraction  angle  of  the  far 
field  propagation  after  the  laser  beam  passes  through  the 
absorption  layers  [6].  When  t>200ms,  with  the  disappearance  of 
the  thermal  lens,  the  divergence  angle  of  the  light  beam  tends  to 
approach  the  diffraction  angle.  This  is  also  consistent  with  our 
theoretical  analysis. 


The  first  draft  of  the  paper  was  received  on  January  25, 
1991;  the  final  revised  draft  was  received  for  publication  on 
April  16,  1991. 
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